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Chapter 2 

PARTIAL DIFFERENTIAL EQUATIONS OF SECOND ORDER 

INTRODUCTION: An equation is said to be of order two, if it involves at least one 
of the differential coefficients r = (d 2 / / d 2 x), s = (d 1 / / dx dy), t =(d 2 z / d 2 y), but 
now of higher order; the quantities p and q may also enter into the equation. Thus the 
general form of a second order Partial differential equation is 

f(x,y,z,p,q,r,s,t) = 0 ...( 1 ) 

The most general linear partial differential equation of order two in two independent 
variables x and y with variable coefficients is of the form 

Rr + Ss + Tt + Pp + Qq + Zz — F ■■■(2) 

where R, S, T, P, Q, Z, F are functions of x and y only and not all R, S, T are zero. 

Ex.l: Solver = 6 x. 


d 2 z 

Sol. The given equation can be written as —7 = 6 x 

Integrating (1) w. r. t. x— = 3x + 0i(y) 

where 0 i(y) is an arbitrary function of y. 

Integrating (2) w. r. t. we get 

xz = x 3 + x 0 i(y) + 02 (y) 
where 0 o(y) is an arbitrary function of y. 

Ex.2. ar — xy 

d 2 z 1 

Sol: Given equation can be written as — = -xy 
Integrating (1) w. r. t., x, we get 



where 0 i(y) is an arbitrary function of y 
Integrating (2) w. r. t., x, 

Z = (a) \ + x 01 ( y ) + 02 ^ y) 


...( 1 ) 

...( 2 ) 


...(1) 

...( 2 ) 
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or z = ^ + x0i(y) + 0 2 (y) 

where 0 2 (y) is an arbitrary function of y. 

Ex.3: Solve r = 2 y 2 
Sol: Try yourself. 

Ex. 4. Solve t = sin(xy) 

d 2 z 

Sol. Given equation can be written as —7 = sin(xy)...(l) 

Integrating ( 1 ) w. r. t., y 

0 = - Q cos(xy) + 0!(x) ... (2) 

Integrating ( 2 ) w. r. t., y 

z = — 07J sin(xy) +y0 2 (x)+ 0 2 (x) 

which is the required solution, 0i, 0 2 being arbitrary functions. 

Exercises :xys — 1 

Sol: We know that s = 

dxdy 

d 2 z 

Therefore xy-= 1 

J dxdy 
d 2 z 1 

or —— - — 

dxdy xy 

Integrating w.r.t., y we have 

dz 1 

- = -log y + f(x) 

Again integrating w.r.t., x we get 

z = log x logy + J fix') dx + F(y) 

Or z = logxlogy + ,g(x) + F(y) 
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Exercises: 2x + 2 y = s 

Sol: The given equation can be written as 

d 2 z 


dxdy 


— 2x + 2y 


Integrating w.r.t., y, we have 


dz 

dx 


= y 2 + 2xy + /(x) 


Integrating w.r.t., x, we have 


z = y 2 x + x 2 y + J /(x) dx + F(y) 
z = y 2 x + x 2 y + p(x) + F(y) 

Exercises:xr + p = 9x 2 y 3 

Sol: The given equation can be written as 

d 2 z 

x—j + p = 9 xy 
dx A 

dp 

=> x — + p - 9x L y A 
dx 

^ 3 JL + l = 9xy 3 

OX X 

which is linear first order differential equation in p 
I. F. is e log x — x 
Multiplying (1) by x we get 


x 


dp p 
Ldx + xJ 


_ g x 2y3 


I 


px = 9 I x 2 y 3 dx 


p* = 9 -y- + /(y) 


px = 3x 3 y 3 + /(y) 


...( 1 ) 
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=> p 


3 x 3 y 3 + f(y) 
x 


dz 

— - 3x 2 y 3 + 
ox 


Integrating with respect to x we get 


/(y) 

x 


z — x 3 y 3 + /(y) log x + F(y) 


Exercises:yt — q — xy 
Sol: Please try yourself. 
Exercises: t — xq = x 2 
Sol: Please try yourself. 
Exercises:r = 2y 2 


Sol: The given equation can be written as 


d 2 z 
dx 2 


2y 2 



Integrating with respect to x we get 


p = 2y 2 x + /(y) 
dz „ 

=> = 2 t x + /(y) 


Integrating we get 

z = y 2 x 2 + J f(y) dx + F(y) 
=> z = y 2 x 2 + x/(y) + F(y) 

Exercises: t = sin(xy) 

Sol: Please try yourself. 

Exercises:logs = x + y 

Sol: The given equation can be written as 
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Integrating w.r.t. x we get 


Integrating w.r.t., y, we get 


dq 

log- = * + y 


dq 

=> -r— = e x e y 
dx 


q = e x e y + fiy) 
dz 

> — = e x e y + /(y) 
dy 


z — e x e y + 


z — e x e y + g(y ) + F(x) 


J /(y) rfy + f(x) 


Exercises:s — t — — 


Sol: Please try yourself. 

Exercises: t + s + q — 0 

Sol: The given equation can be written as 


dq dp dz 
dy dy dy 


Integrating with respect to y, we get 


q + p + z — f{x) 

=> p + q = fix) — z 


It is of the form Pp + Qq = R 
Its auxiliary system is 

dx dy dz 

1 1 f(x)—z 

From first two fractions of (1) we get 
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dx — dy 


=> x — y — a 

From first and third fractions of (1) we get 


dx 


dz 


1 f{x) - z 
> [fix) — z]dx = dz 
dz 

*Tx = f{x) ~ z 

dz 


It is first order linear differential equation in z 
Its integrating factor is e^ dx —e x 
Therefore ze x = / fix) e x dx 


ze 


— J f(x) e x dx + f(y) 


Exercise:! + s + q — 1 
Sol: Please try yourself. 

Exercise: Find the surface passing through the parabolas, 

y 2 = 4 ax, z — 0 
and y 2 = —4ax, z = 1 

and satisfying the equation xr + 2p = 0. 

Sol: The given second order partial differential equation is 

xr + 2p — 0 


+ -p=0 

OX X 


It is first order linear differential equation in p. 

Its integrating factor is e^ dx = e 2logx = e log * 2 = x 2 


...( 1 ) 
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From (1) we get 


x2p = I 0dx 

*V-/od* + /M 


/GO 

dis¬ 


integrating w. r. t. x we have 

z = -r/GO + F(y) 


Using the given condition z = 0, x — —, in equation (2), we have 

1 4 a 

4 CL 

0 = /(y) +F(y) 

y£ 


or F(y) 


4a/(y) 


Also for z = 1, and x —-, we have from (2) we have 

4 a 

4 CL 

i = —/GO + F(y) 
y 


Using (3) we get 


or 


_ 4a/(y) 4a/ (y) 

-L 9 I 9 


=> 1 = 


8a/(y) 

y 2 


Substituting /(y), in (3) 


y 

/w =s 


4a y 2 


F W = 2 


- ( 2 ) 


...(3) 
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Therefore from (1) we get 


_ -y 2 1 

8 ax 2 


Which is the required surface passing through the parabolas. 

Exercise: Find the surface satisfying t = 6x 3 y containing the two lines 
y = 0 = zand y = 1 = z 
Sol: The given 2 nd order PDE is 

t — 6 x 3 y 



Integrating w. r. t., y, we have 


R = 


6x 3 y 2 

2 


+ /(x) 


dz , „ 

— - 3x 3 y 2 + f{x) 
oy 


Integrating w. r. t., y, 

3x 3 y 3 

z = —-— + y/(x) + F{x) 
=> z — x 3 y 3 + y/(x) + F(x) 

Using given conditions y = 0 = z, in (1), we have 

0 = 0 + 0 + Fix) 

=> F(x) = 0 

Also using y = 1 = z in equation (1) we get, 

1 = x 3 + fix) + Fix) 
Using (2), we get 1 = x 3 + fix) + 0 
fix) — 1 — x 3 

Using (2) and (3) in (1) we get 


...( 1 ) 


• ••( 2 ) 


• ••(3) 
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z — x 3 y 3 + y( 1 — x 3 ) 

Which is the required surface containing the two lines. 

Exercise: Find the surface satisfying r + s — 0, and touching the elliptic paraboloid 
z = 4x 2 + y 2 along the surface of plane y = 2x + 1. 

Sol: From the given equation we have ^ ^ = 0. 


Integrating with respect to x, we have 

V + q = f(y ) 

Now, the auxiliary system is 

dx _ dy _ dz , 1 . 

Taking first two fractions we get 

dx dy 

T“T 


Integrating we get 


x — y + a 

=> x — y — a •••(2) 

Also from 2 nd and 3 rd fractions of (1), we get 

dy dz 

1 /GO 

=> dz = /(y)dy 
=> z = (y) + b 


or z = <p(y) + F(a) 


=^> z = <p(y) + F(x — y) 

• ••(3) 

From (3), we get 


v = j; = F \ x -y') 

...(4) 

q =^ = <p' Cy) y) 

...(5) 
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Since z = 4x 2 + y 2 


dz 

••• p — — = 8x 

“ dx 

• ••(6) 

o dz 

& q=37 = 2 y 

• ••(7) 

From (4) and (6) 


1 

II 

00 

• ••(8) 

From (5) and (7) 


<p' (y) ~P\x~y) = 2y 

• ••(9) 


Adding (8) and (9) we get 

(p (y) = 8x + 2y 
8 

= 2 (y -!) + 2 y 
= 6y — 4 

Integrating w. r. t., y, we get 

<p(y) = 3y 2 — 4y + b ...(10) 

Also, from (8) 

—F' (x — y) = 8x = — 8(y — x — 1) = 8(x — y + 1) 

Integrating w. r. t., (x — y) we get 

—F(x — y) = 4(x — y) 2 + 8(x — y) + c •••(11) 

Substituting (10) and (11) in (3) we get 

z = 3y 2 — 4 y + b — 4(x — y) 2 — 8(x — y) + c 
= —4x 2 — y 2 + 4y — 8x + 8xy + d 
From the given condition, 

4x 2 + (2x + l) 2 = —4x 2 — (2x + l) 2 + 4(2x + 1) — 8x + 8x(2x + 1) + d 
=> 8x 2 + 2(2x + l) 2 = 4(2x + 1) — 8x + 8x(2x + 1) + d 
=> 8x 2 + 8x 2 + 2 + 8x = 8x + 4 — 8x + 16x 2 + 8 x + d 
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=> d — —2 


Therefore z = —Ax 2 — y 2 + Ay — 8x + 8 xy — 2 

which is required surface. 

Exercise: Show that the surface satisfying r = 6x + 2 and touching z — x 3 + y 3 
along its section by the plane x + y + 1 = 0 is z = x 3 + y 3 + (x + y + l) 2 . 

Sol: Try yourself. 


Partial differential equations with constant coefficients: 

We know that the general form of a linear partial differential equation 


d n z 


+ Ar,n 


d n z 


d n z 

71 dx n n ~ 1 dx n ~ 1 dy " T " n ~ 2 dx n ~ 2 dy 2 


+ "' + y4l P = ^ (X ' y) 


...( 1 ) 


Where the coefficients A n , A„_ 1( A n _ 2 , ... , A 1 are constants or functions 

of x and y. If A n , A n _ lt A n _ 2 , ... , A 1 are all constants, then (1) is called 

a linear partial differential equation with constant coefficients. 

We denote and ^ by D (or D x ) and D (or D y ) respectively. 


Therefore (1) can be written as 

[A n D n + A n _ 1 D n ~ 1 D' + A n _ 2 D n ~ 2 D' 2 + ... + A x D’ n ]z = /(x,y) ... (2) 

or (p{D,D')z = /(x,y) 


The complementary function of (2) is given by 

\A n D n + A^D^D' + A n _ 2 D n ~ 2 D' 2 + ••• + AiD' n ]z = 0 ...(3) 

or (p(D, D')z — 0 

Let z = F(y + mx) be the part of the solution 

Dz — ^ — mF' (y + mx) 

D 2 z = = m 2 F" (y + mx) 
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And 


D n z - = m n F n (y + mx) 


D z — — F (y + mx) 

D 2 z — = F" (y + mx) 


D n z — = F n (y + mx) 

Substitute these values in (3), we get 

[A n m n + An^m 71-1 + A n _ 2 m n ~ 2 + . . . + AjF^Cy + mx) = 0 

which is true if 'm' is a root of the equation 

If m t , m 2 , m n , are distinct roots, then complementary functions is 
z = <pi( y + mi x) + cp 2 (y + m 2 x) + .. . +cp n (y + m n x) 

where (pi, (p 2 , ■ . (p n are arbitrary functions. 

••• (p(D,D')z = 0 

we replace D by m and D' by 1 to get the auxiliary equation from which we get 
roots. 


Linear partial differential equations with constant coefficients 


Homogenous and Non homogenous linear equations with constant coefficients: A 

partial differential equation in which the dependent variable and its derivatives appear 
only in the first degree and are not multiplied together, their coefficients being 
constants or functions of x and y, is known as a linear partial differential equation. 

The general form of such an equation is 


z 

A °^ + A 


d n z 


1 a„a ~n 1 1 Qy2g x 


B i 


dy 2 dx n 3 


dyd x r 

+ . . . 


d n z d n zl d n ~^z d n ~^z 

-dx n ~ 2 n dy n . ® dx n ~ 1 ^ dyd x n ~ 2 


+ B„ 


dy n 


t] + [ 


. , dz , „ , dz 

m «57 + m ^ 


+ z = f(x,y) 


dydx n 2 

...( 1 ) 
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where the coefficients A 0 , Ai, . . . A n , B 0 , Bi,. . ., B n i, M 0) Mi and N 0 are all 
constants, then (1) is called a linear partial differential equation with constant 
coefficients. 

For convenience and will be denoted by D and D' respectively. 


2.1. Short method for finding the P.I. in certain cases of F(D,D) z = f(x, y) 

2.1. a. Short method I. when f(x, y) is of the form 0 ( ax + by ) 

Ex.l. Solve (D 2 +3DD' +2 D' 2 )z = x + y 

Sol. The Auxiliary equation of the given equation is 
m 2 +3m +2 = 0 giving m = -1 ,-2 

therefore C.F. = 0i(y-x) + 0o(y-2x), 0i, 02 being arbitrary functions 


Now 


P.I. 


l 


D2 +3DD' +2D'2 
1 


(x + y) 


1 2 +3.1.1+2.1 2 


// vdvdv, where v= x+y 


fO 


lr 1 , .3 

-= — (x + y) 

6 6 36 v : ’ 


Hence the required general solution is z = C.F. + P.I. 


or z = 0i(y-x) + 0 2 (y-2x) + — (x + yf 


Ex. 2. Solve (2 D 2 -5 DD' + 2D' 2 ) z = 24 (y - x ) 
Sol: Try yourself 

Ex.3. Solve (D 2 +3 DD' +2D' 2 )z = 2x + 3 y 
Sol: Try yourself 

2.1.b.Short method II. 

When f(x,y) is of the form x m y n or a rational integral 
Ex.l. Solve (D 2 -a 2 D ,2 )z = x 

Sol. Here auxiliary equation is nT -a =0 so that m = a, -a 
Therefore C.F. = 0 1 (y + ax) + 0 2 {y — ax), 


...( 1 ) 
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0 1; 0 2 being arbitrary functions. 

i i 


Now P.I. = 


D 2 -a 2 D' 2 


(x) = 


D 2 1 — 


-(£)] 


(x) 


1 

D2 


“ 2 ©] 


= X[1 + a’- (D -/D 2 ) +...]x 


1 


D2 

„3 


Hence the required solution is z = C.F. + P.I. 


Z = 0i(y+ a x ) + 02 (y-a x) + — 

o 


Exercise: 1.2r + 5s + 2t = 0 


Sol: It is a second order pole with constant coefficients, we have 

d 2 z d 2 z d 2 z 
^ dx 2 + ^ dxdy + ^ dy 2 

or (2D 2 + SDD’ + 2D’ 2 )z = 0 

Now the auxiliary equations is given by 

(2m 2 + 5m + 2) = 0 

=> m = — - , -2 

Therefore the complementary function is z — (p 1 (y — ^-x^ + (p 2 (y ~ 2x) 
which is required solution. 

Exercise: 2.r = a 2 t 


Sol: Try Yourself (Ans: z — cp t (y + ax ) + cp 2 (y — ax) ) 


„ . _ d 3 z „ a 3 z 

Exercise: 3.—— — 3 - 

dx 3 dx 2 dy 


+ 2 


d 3 z 

dxd 2 y 


Sol: It is a third order pole with constant coefficients, we have 


...( 2 ) 
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d 3 z d 3 z d 3 z 

—- 3— -h 2-■— = 

dx 3 dx 2 dy dxd 2 y 

or (D 3 — 3D 2 D’ + 2DD' 2 )z — 0 

Now the auxiliary equations is given by 

m 3 — 3 m 2 + 2m = 0 
=> m(m 2 — 3m + 2) = 0 
=> m(m — l)(m — 2) = 0 


=> m = 0, 1 , 2 

Therefore the complementary function is z = (p 1 (y) + q) 2 (y + x) + (p 3 (y + 2x) 
Which is required solution. 

_ • a 9 3 z . d 3 z . . d 3 z , d 3 z - 

Exercise: 4 .—7 — 6 —+11 ——5 - 6 —r — 0 

dx 3 dx 2 dy dxd 2 y dy 3 

Sol: Try yourself 

Exercise: 5. 25r — 40s + 16t = 0 


Sol: Try yourself 

Exercise: 6.(D 4 — D' 4 )z = 0 

Sol: The auxiliary equation is given by 

m 4 — 1 — 0 


=> ( m 2 — 1 )(m 2 + 1 ) = 0 
=> (m — 1 )(m + 1) (m — i)(m + i) = 0 
=> m — 1, — 1, i, —i 


Therefore the complementary function is 

z = ^(y + x) + <p 2 (y -x) + <p 3 (y + ix) + <p 4 (y - ix) 
which is required solution. 

Exercise: 7.(D 3 — 4D 2 D' + 4 DD’ 2 )z — 0 
Sol: Try yourself. 
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Exercise: 8 .(D 2 — 2DD' + D' 2 )z = 12 xy 

Sol: The auxiliary equations corresponding to these linear system of equations is 
given by 

m 2 — 2m + 1 = 0 


=> (m — 1) (m — 1) = 0 

=> m — 1, 1 


Therefore the complementary function is 

2 = fi(y + x) + xf 2 (y + x) 


Also, P. I. 


(D-D') 2 


12xy 



xy 


_ 12 
_ ~D 

= 12 



= 2x 3 y+ x 4 


Therefore the complete solution is z=C. F. + P. I. 
i.e., z — f t (y + x) + xf 2 (y + x) + 2x 3 y + x 4 
Exercise: 9.(2 D 2 — SDD' + 2 D' 2 )z — 24(x — y) 
Sol: Try yourself. 

Exercise: 10.(D 3 — D' 3 )z — x 3 y 3 
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Sol: The auxiliary equations is 


m 3 — 1 = 0 

(m — 1 )(m 2 + m + 1) = 0 

-1 + iV3 — 1 — iV3 


m = 1, 


Complementary function is z = (p t (y + x) + (p 2 (y + + <jt >3 (y + 

-l-iV3 


■ X ) 


P. I. 


3,.3 


—rrx y 

D 3 —D 3 J 


3,.3 


D 3 1- 


-©] 


TfX y 


1 

£> 3 

1 

D 3 


1 - (-)’] 1 x 3 y 3 




3 3 

x y 


1 3 3 - \ 3 3 1 fD'\ 3 3 1 

= ^3 *' y + (-) x 'y + (-) x 'y + 


V 2 


1 

D 3 

1 

D 3 


^[x 3 y 3 +^3x 3 y 2 + 0+ . . . ] 
^3 [x 3 y 3 + p- 6 x 3 y] 

x 3 y 3 + ^- 6 x 3 J 
3 3 , 1 ,-* 4 l 

[x y + ^76—J 

1 f 3 3 . 1 r x5 l 
= —7 x y + -6 — 

D 3 L * D 20j 

= p[*V + f^] 

J_ | x 4 y 3 jcM 
fl 2 i 4 ^ 140 J 

= I f^v + je!_1 

d L 20 1120 J 
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r x 6 y 3 x 9 


120 10080J 

The complete solution is C.F. + P.I. 


2 = <Pi(y + *) +(p2 y + 


-1 + iV3 


■x) + (p 3 y + 


-i - ;V3 


■x + 


x 6 y 3 


+ 


X' 


120 10080 


Exercise: Find the real function V of x and y reducing to zero when y — 0 and 
satisfying — + — = -4tt(x / + y z ) 

Sol: We have to find the P. I. only 


P. I. 


d 2 +d' 2 

i 


— [—4tt(x z + y 2 )] 


D 2 


(»S) 


77 t[- 47 t(x 2 + y 2 )] 


—An 

D 2 

1 + 

D ' 2 1 1 
~D 2 \ 

(x 2 + y 2 ) 


II 

to I 

1 - 

— + 1 

D 2 ' 

® 2 + --] 

(x 2 + y 2 ) 

—An 

D 2 

(x 2 

+ y 2 ) 

+ y 2 ) + 0 ] 

—An 

D 2 

(x 2 

+ y 2 ) 



II 

to 1 X 

"(x 2 

+ y 2 ) 

-F®] 


II 

to I 

M | si 

"(x 2 

+ y 2 ) 

-j(2x)] 


—An I 

D 2 ' 

[(x 2 

+ y 2 ) 

- (x 2 )] 


—An , 

D ' 

\xy 2 

] 




-47T 


r x 2 y 2 1 


2 . 
2,.2 


= —2nxy 

Theorem: If it 1( it 2 , ii 3 , . . . ,u n are the solutions of the homogeneous linear 
PDE F{D, D')z — 0, then Xr=i c r u r where c r r s are arbitrary constants, is also a 
solution. 
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Proof: Since u r , r=l, 2, 3, . . n are solutions of the PDE F(D,D')z — 0 
So, we have u r is one of the solutions i.e., 

F(D,D r )u r = 0 , r = 1, 2, . . . , n 

••• F(D,D )(c r it r ) = c r 'F(D, D )(u r ) 
and F(D,Z/)(5> r ) =XF(D,D')(u r ) 
for any set of functions u r , we have 


( n \ n 

V c r u r ) = / F(D,D )c r u 

hi ) hi 


n 

^ c r F(D,D )u r 

r=1 


= o 


Therefore Y,r=i c r u r acts as a solution for the homogeneous system. 


Reducible and irreducible: 

If an operator F{D, D') can be expressed as a product of linear factors, it is said to be 
reducible. If it can not be factorised, then it is said to be irreducible. 

Theorem: If a r D + /?,.D + y r is a factor of F(D, D') and (p r ), then 

u r = exp. (— —) (p r ([3 r x — a r y ) for a r is a solution of the equation F(D,D')z — 
0 . 


Proof: The given equation is F(D,D')z — 0 ... (1) 

In order to prove u r = exp. (— —) cp r (/? r x — a r y) ; a r ^ 0 . . 

is a solution of (1), we have to prove F(D, D')u r = 0 
Diff. eq.(2) w.r.t. x and y, we get 

Y r ( y r x\ 

Du r - - u r + P r exp. - (p (B r x — a r y ) 

a r V a r / 


( 2 ) 
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And 


D'u r — —a r exp. j <p'(/? r x — a r y ) 

, / y r x\ , 

••• ( a r D + /?,. D + Yr)u r — —y r u r + a r (3 r exp. ^ - J (p (/? r x — a r y) 


—a r /3 r exp. ( — —) (p (/? r x — a r y ) + y r u r — 0 . . . (3) 

Since ( a r D + /3 r D + y r ) is a factor of F(D, D') 

Therefore F(D, D')z — g(D, D')(a r D + /3 r D + y r )z, using (3), we get 

F(D,D')u r — 0 

Therefore u r is a solution of F(D, D')z — 0 


Solution of Reducible Equations: 

Let F(D, D')z = f (x, y) ... (1) 

be a partial differential equation. Since (1) is reducible therefore 

n 

F(D, D')z — J”|(a r D + /3 r D + y r )z 

r=1 

If z satisfies ( a r D + jS r D + y r )z = 0, r = 0,1, 2, . . . , n, then it gives us 
complementary function 

Now a r ~r + /3 r ~r + y r z — 0 

ox oy 

The subsidiary system is 


dx dy dz 
a r /? r y r z 


From the first two members 


/3 r x — a r y — c r 

From first and last members we get 


dz y r 

— =- dx 

z a r 
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= log Pr exp [ eA = fir] 

=> z = (r) exp — xj 

z - <p(/? r x — a r y) exp — xj 


dz y r 

T = ~K dy 

=> z = <p(/? r x) exp —yj 

Example: Let ( a r D + (3 r D' + y r )z 1 = 0 where z x = (a r D + (3 r D + y r )z 

z 1 — <p(/? r x — a r y) exp — xj 

=> ( a r D + p r D + y r )z = <p r (/? r x — a r y ) exp — xj 

dz dz ( y r \ 

=> a r — + /?,- — = (p r (J] r x - a r y ) exp-x ) — y r z 

dx dy \ a r / 


Auxiliary system is 


dx dy dz 

a r fir cp r (/? r x — a r y) exp ( — — Yr z 


From first two we get 


dx dy 
a r /?,. 

=^> /? r dx = a r dy 
/?,. x — a r dy — c r 


From first and third we get 
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dx dz 

a r (p r (/?,. x - a r y) exp ( — x) - y r z 

dz (p r (/?,- x - a r y) exp (- ^ *) - V r z 
dx a r 


dz y r z <Pr(firX-a r y)ex p(-^x) 


dx + a 


a r 


Yr x 


Here I. F. is e a r therefore the above equation can be written as 

YzL\ (p r (p r x - a r y) 


( LL-—\ 

Ize «r J = 


a, 


YrX 1 

ze a r 


1 f 

— — (p r {f! r x — a r y)dx + /?,. 
cc r 


y r X 


y r X 


=$ ze a r = e a r [cp r Q3 r x — a r y) + t p r Q3 r x — a r y )] 
Example: If z = e ax+by 


Then F(D,D')z = F(a,b)e ax+by 


z acts as the solution of F(D, D r )z, where F{D, D')z is reducible if F(a, b ) = 0. 


„ . d 3 z d 3 z 

Exercise:— 7 — 2 — 7 — 

dx 3 dx 2 dy 


d 3z , 2 — = e x+y 

dxdy 2 dy 3 


Sol: The given differential equations can be written as 

(D 3 - 2D 2 D - DD’ 2 + 2D' 3 )z = e x+y 
Auxiliary equations m 3 — 2 m 2 — m + 2 — 0 

=> (m — l)(m 2 — m — 2) = 0 
=> m — 1, — 1, 2 

Therefore the C. F. is 

2 = A(y + x) + fi(y-x) + A(y + 2x) 

p t = _ - _ px+y 

D 3 —2D 2 D’ —DD' 2 +2D' 3 
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- _1_ e x+y 

D 2 (D-2D')-D ,2 (D-D') 

- _ l _ e x+y 

(d 2 —D' Z )(D—2D') 


- _ l _ e x+y 

(D-D')(D+D')(D-2D') 


_ _ l _ e x+y 

(D-D) (1+1) (1—2) 


- 1 c x+y 

-2 (£>-£>') 

Now let w — ; e x+y 

D-D' 

=» ( D-D')w = e x+y 

The auxiliary equations are 

dx dy dw 
e x+y 

From first two members we have 

dx dy 

=> dx + dy = 0 
=> x + y — c 

From first and third member we get 

dx dw 

T ~ e x+y 

dw 

=> dx = — 
e c 

=> dw — e c dx 


=> w = e c x 


=> w = xe x+y 


Therefore the particular integral = —j — — ^xe x+y 
Hence the complete solution is z — C.F + P.I 
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z — fi(y + x) + fi(y — x) + /i(y + 2x) H— -xe x+y 


Laplace Equation 


V 2 z = 


a 2 a 2 \ 

dx 2 + dy 2 ) Z 


Exercise: Find the solution of the equation V 2 z = e -x cos y 
Which tends to 0 as x oo and cosy for x = 0. 

Sol: The given pde is V 2 z = e _x cos y 

a 2 z a 2 z 

^3^ + ^ =e " C0Sy 


(D 2 + D' 2 )z = e x cosy 


On comparing with F{D, D') — D 2 + D' 2 and /(x, y) = e x cos y 
Let z = e ax+by be the solution of (1) 

(D 2 + D' 2 )z = a 2 e ax+by + b 2 e ax+by 
= ( a 2 + b 2 )e ax+by 

Where a 2 + b 2 — 0 — F(a, b ) 

Therefore the complementary function is C. F. — A r e ax+hy 

A r 's being the constants and a 2 + b 2 = 0 

1 _ 

Also P. I. = —r—— e x cos y 

d 2 +d' 2 j 


cos y— i e 


x cos y — e" 

J 2D 


■cosy e" 


Therefore the complete solution is 


z = ^ A r e ax+by — - cosy e" 
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Using z 0 as x ^ oo, we write 
a r — —A r where A r > 0 
Since a 2 + b 2 — 0 

=> b r — ±iy[af 

— ~f~ i A < ■ 


z 


00 

-\cosye~* 

r =0 


00 

= ^ B r e~ XrX cos(A r y + e r ) — - cosy e 

r=0 

Using the boundary condition 


cosy = 


00 

B r cos {A r y + e r ) 

r =0 


Where B r — 1 and A r — 1 for r — 0 and B r — 0 and A r — 0 
Therefore z = cos y e -x — - cos y e~ x is the required solution. 

Exercise: Show that the equation —4 + 2k — — c 2 —-r 

M at 2 at ax 2 

Possesses solution of the form Y. C r e~ kt cos(a r x + c r ) cos(w r £ + S r ) 
Where C r , a r , e r , a) r , 5 r are constants and co r = S 2 c 2 — k 2 

Sol: Let ^ = D; ^ = D'\ 

Therefore (D 2 + 2 kD)y — c 2 D' 2 y 
or (D 2 + 2/cD - c 2 D' 2 )y = 0 

It is irreducible, therefore let y = e at+bx 
=> Dy — ae at+bx and D 2 y — a 2 e at+bx 
Similarly we get 

D'y — be at+bx and D' 2 y — b 2 e at+bx 


for r =£ 0 


...( 1 ) 
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Therefore form (1) we have 

a 2 e at+bx + 2 kae at+bx + c 2 b 2 e at+bx = 0 

=> (a 2 + 2 ka + c 2 b 2 )e at+bx — 0 

=> a 2 + 2 ka + c 2 b 2 = 0 

—2k + V4/c 2 + 4 b 2 c 2 

=> a =- 

2 

=> a — —k ± yjk 2 + b 2 c 2 
In general a r — — k + yjk 2 + b 2 c 2 

If b 2 — —a 2 
Then a r — —k + J k 2 — a 2 c 2 
= —k + im r 

Where a) 2 — a 2 c 2 — k 2 


Therefore y — e b,x e a,t 

_ g—xtg—ia) r tg—ia r x 
00 

y = I c r e xt cos(a r x + e r ) cos(co r t + 8 r ) 

r =0 

Where C r , a r , e r , a> r , 8 r are constants and n> r = S 2 c 2 — k 2 

Exercise:l Ifz = /(x 2 — y) + g(x 2 + y), where / and g are arbitrary constants, 
prove that 

d 2 z 1 dz d 2 z 

-— 4 X 2 - 

dx 2 xdx dy 2 

Exercise:2 Find the solution of 

d 2 z d 2 z 
d^ 2 ~d^ 2=x ~ y 
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Sol: The given equation is 

(D 2 — D 2 )z — x — y 

The auxiliary equation is m 2 — 1 = 0 so that m = +1 
Therefore C. F. is A (y + x) + f 2 (y — x) 

P. I. is - 1 , 2 (x - y) 

D 2 —D 



= -y) + 0 + 0 ■ ■ ■ ] 

__ i 

_ D 


3 

Therefore the complete solution is z = C. F. +P. I. = A (y + x ) + A (y — x ) H- 

? 

2 

Exercise:3 Find the solution of 

d 4 z d 4 z d 4 z 

dx 4 + dy 4 ^ dx 2 dy 2 

Sol: Please try Yourself 
Exercise:4 Show that the equation 


x 

~2 


yx 


x 

T 


yx 
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d 2 z 1 dz 
dx 2 k dt 

Possesses solution of the form 


00 

r =0 


Sol: Please try Yourself 

Exercise:5 Find the P. I. of the following PDE’s 

(a) ( D 2 — D')z — 2 y — x 2 

(b) (D 2 - D')z = e 2x -y 

(c) r + s — 2t — e x+y 

(d) r — s + 2q — z — x 2 y 2 


d 2 z 


d 2 z 

dx 2 dxdy 


. d 2 z 

~ b Hy 2 


y cosx 


Sol: (a) P. I. is ^^(2 y - x 2 ) 


1 

D 2 


D 2 

1 

D 2 


1 -^ 
1 D 2 


(2 y — x 2 ) 


1 - 


D' 


n _1 


D 2 


D' D' 2 

1 + W + ~D* + ' ' ' 


(2 y — x 2 ) 

(2 y — x 2 ) 


1 

D 2 


D' D' 2 

(2 y - * 2 ) + ^2 ( 2 y _ + - x2 ) + ■ 


1 

D 2 


(2y — x 2 ) + ^(2) 


^2 (2y-x 2 ) +^(2x) 


D 2 


[(2y — x 2 ) + x 2 ] 
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= >y] 

= [2xy ] 

2 

= x z y 


Which is required particular intergral 

(b) Please try yourself. 

(c) Please try yourself. 

(d) P. I. is —;- r -t—( x 2 y 2 ) 

v ’ D 2 -DD +2 D -1 v J 


D 2 

[l_^ + 

2D' 

1 ' 

1 D + 

D 2 

d 2 _ 


(x 2 y 2 ) 


1 

1 _ 


2D' 

1 ' 

_|_ 

D 2 

_L 

u 

D 2 

D 2 


-l 


(x 2 y 2 ) 


1 

D 2 



2D' 

1 ' 

_|_ 

\D 

D 2 

D 2 



2D' 

1\ 2 


u 

D 2 

D 2 ) 



(x 2 y 2 ) 


l 

D 2 


[(x 2 y 2 ) + i (2x 2 y) - ^ (2x 2 y) + ^ (x 2 y 2 ) + ^ (x 2 y 2 ) + 

4^(x 2 y 2 )+^(x 2 y 2 )-4^(x 2 y 2 )-^(x 2 y 2 )+^(x 2 y 2 )+ . . 




= p [(x 2 y 2 ) + ^ (x 3 y) - ^ (x 4 y) + i (x 4 ) + ^ (x 6 ) - ^ (x 5 ) 
^(x 6 y) +^(x 5 y)] 


x^ >x^y x 4 y 2 x 6 y >x^y x 2 y 

2240 315 + 180 + 30 + 12 120 2570 630 


which is the required particular integral. 
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Chapter-3 

Classification of second order partial differential equation 

Definition: A second order partial differential equation which is linear w. r. t., the 
second order partial derivatives i.e. r,s and t is said to be a quasi linear PDE of 
second order. For example the equation 

Rr +Ss +Tt + f(x,y,z,p,q) — 0 ■■■■(!) 

where f(x, y, z, p, q) need not be linear, is a quasi linear partial differential equation. 
Here the coefficients R, S, T may be functions of x and y, however for the sake of 
simplicity we assume them to be constants. 

The equation (1) is said to be 

(i) Elliptic if S 2 -4RT < 0 

(ii) Parabolic if S 2 -4RT =0 

(iii) Hyperbolic if S 2 -4RT >0 

BOUNDARY VALUE PROBLEMS: The function v in addition to satisfying the 
Laplace and Poisson equations in bounded region R in three dimensional space, 
should also satisfy certain boundary conditions on the boundary C of this region. 
Such problems are referred to as Boundary value problems for Laplace and poison 
equations. If a function f G C n , then all its derivatives of order n are continuous. If f G 
C°,then we mean that f is continuous. 

There are mainly three types of boundary value problems for Laplace equation. If f 
G C 0, and is prescribed on the boundary C of some finite region R, the problem of 
determining a function 0 (x, y, z) such thatV“0 = 0 within R and satisfying 0 = f on 
C,is called the boundary value problem of first kind or Dirichlet problem. The second 
type of boundary value problem (BVP) is to determine the function 0 (x, y, z) so that 

V 2 0 = 0 within R while — is sepecified at every point of C, where —is the normal 

on on 

derivative of 0. This problem is called the Neumann problem. 

The third type of boundary value problem is concerned with the determination of the 
function 0 ( x,y,z ) such that V“0 = 0 within R,while a boundary condition of the 

form h 0 = f,where h> 0 is specified at every point of the boundary C.This is 
called a mixed boundary value problem or Churchill’s problem. 
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Separation of variables method: 

The method of separation of variables is applicable to a large number of classical 
linear homogenous Equations. The choice of the coordinate system in general 
depends on the shape of the body. Consider a two dimensional Laplace equation 


V 2 0 

S 2 (u) d 2 (u) 

dx 2 dy 2 

-(1) 

We assume that 

u(x,y ) = X(x) Y (y) 

...(2) 

Equation (1) and (2) provide us 


X - = Y - = k 

X Y 

(separation parameter) 



Three cases arise : 

Case I: Let k>0.then k = p“, p is real we get 
- p 2 X = 0 and ^ - p 2 X = 0 

ax z ay z 

which imply that X = Ci e px + Co e' px 
and Y = C 3 cos py + C 4 sin py 
then solution is 

u(x,y) = (Ci e px + C 2 e' px ) (C 3 cos py + C 4 sin py ) ...(3) 

Case II: let k = 0 then y-y — 0 and = 0 

ax z ay 

Which provide us X = C 5 X + C 6 and Y = C 7 y + Cg 

The solution is therefore u(x,y) = (C 5 X + C (,)( C 7 y + Cg) ...(4) 

CASE III: let k < 0 then k = -p proceeding as in case I, we obtain 

u(x,y) = (C 9 COS px + C 10 sin px) (C 11 e py + C 12 e' py ) ...(5) 

In all these cases Cj ( i = 1, 2, 3, . . . ,12) are integration constants,which are 
calculated by using the boundary conditions. Lor example, consider the boundary 
condition 

U(x,0) = 0, u(x, a) =0, u(x, y) -» 0, as x -> 00 
Where x > 0 and 0< y < a. 
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The appropriate solution for u(x, y) by the methods of separation of variables 
obtained above in this case is 

U(x, y) = (Ci e px + C 2 e' px ) (C 3 cos py + C 4 sin py ) ...(6) 

Since u(x,y) -* 0, as x -* oo, we have 
Ci = 0 V y 

(x,y) = C 2 e' px (C 3 cos py + C 4 sin py ) 

As u(x,0) = 0, we get 

C 2 e' px C 3 = 0 => C 3 = 0 [ because C 2 A 0A e~ px ’ V x ] 

U(x, y) = A e' px sin py, A = C 2 C 4 

Now u(x, a) = 0 => A e' px sin pa = 0 => sin pa = 0 [v A A 0 ] 

=> pa = n 7 r,nG I => p = n 7 r / a, n = 0, + 1,... 
u(x, y) = X A n e‘V x/a sin 

A n being new constant. 

This is the required solution in this case. 

Ex. Show that the two dimensional Laplace equation VC V = 0, in the plane polar 
coordinates r and 9 has the solution of the form (Ar 11 + Br 11 ) e (+ in0), 

where A and B are n constants. Determine V if it satisfies 

Vi 2 V = 0 in the region 0< r < a, 0< 9 <2n and 

(i) V remains finite as r 0 

(ii) V = Hn Cn cos(n0), on r = a. 

Sol: Try yourself. 

Laplace equation in cylindrical coordinates: 

^ 2 t/ d 2 V 1 dV 1 d 2 V d 2 V . 

V ^2 r dr r 2 dQ 2 q z 2 

Let V = R(r)0(0)Z(z) be the solution 
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dV dR A d 2 V ^d 2 R 

— = 0Z — and —— = 0Z —— 
dr dr dr' dr' 


d 2 v _ 07 d 2 & 
do 2 do 2 ' 


<8 ^ 

dz 2 dz 1 


■ v *v=ez^? + iez^ + -LjE«!? + i»«:f 

v dr 2 r dr r 2 d# 2 dr 


^ 1 1 dR 1 J 2 0 1 </ 2 Z „ 

Or- - +-+ —- T +- T = 0 

R dr 2 r dr r'& dO 2 Z dz 2 


T 1 d 2 Z 2 , d 2 Z 2 _ _ 

Let ~— = m , then —~ m Z = 0 

Z dz dz 


:.Z = , 


r 1 d 2 0 

Now let-- 

0 dO 2 


Cl-O 2_ n 
— T + n 0 = 0 
dO 2 n 




d~R 1 dR 2 n 

Now — + - — + R m -j =0 

dr r dr „ 

v r J 


which is Bessel’s equation 


its solution can be written as 


R(r)= A,J ,hn) + B,Y ,hn) 

Therefore V(r,0,z) = { A J (m, )+ B,Y.ittt >/" Y* 


Home Assignments 


Exercise: Solve the PDE 


8 2 (j> 1 d<f> 1 d 2 (j) 

-Yr +—- + ——T = ° 
dr' r dr r~ d0~ 


Subject to the condition v = — = 0 at r=a 

dr 
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And v r = U cos0, 


v e 


1/M 

/ r dd 


-U sin# as r-» oo 


Exercise: Solve the BVP 


V 2 u = 00 < r < 10 
0 < d < n 


Subject to the conditions 

400 

u(lO,0) =-( nd-d 2 ) 

n 

u(r, 0) = 0 = u(r,n ) 

And it(0, 0) is finite 

Exercise: Show that the solution of the Cauchy problem for the Laplace equation 

d 2 u d 2 u 
dx 2 + dx 2 

Subject to the condition u(x, 0) = 0 

1 

u v (x, 0) = — sin nx 
7 n 


Where n is a positive integer and 

1 

it(x,v) = r sinh ny sin nx 
n z 

Interior Dirichlet problem for a circle 

The Dirichlet problem for a circle is defined as follows: 

To find the value of u at any point in the interior of the circle r = a in terms of its 
values on the boundary such that u is the single valued and continuous function 
within and on the circular region and satisfies the equation V“u = 0 ; 0< r < a subject 
to u(a, 0 ) = f( 0 ); 0 < 6 < 2n 
We have 

V 2 u = d 2 u / dr 2 + 1/r — + l/r 2 5 2 u ldd 2 = 0 

or 

We know that 

u(r,0) = Xn=o(/nr n + D n r" n ) (A n cos nd + B n sin nd ) 
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Since the function u is defined for all the values within and on the boundary of a 
circle 

Therefore for r = 0 u(r,0 ) exists only for D n = 0 V n 
Thus u(r,0 ) = H?T=o T a (A n cos n0 + B„ sin n0 ) 

= a 0 / 2 + Y,n= l r " (A n cos n0 + B n sin n0 ) where A 0 = a 0 / 2 
= ao / 2 + Y.n=l r " (a n cos n0 + b n sin n0 ) 

Put u(a,0 ) = f(0) 

Therefore f(0) = = ao / 2 + S“=i a " (a n cos n0 + b n sin nd ) 

This is the full range Fourier series (i.e. a n , b n A 0 ) 

Now a 0 = 1/7T f Q f(cp)dcp 

2.71 

a n = l/7ra n f f(cp ) cos rnp dcp 

7tt 

b n = l/7ra n f Q f(cp) sin ncp dcp 

2 77 2 77 

u(r,0 )-l/2n / f((p)dcp + l/n'Zn=i(j/ a Tf 0 cos nd cos ncp + sin nd sin ncpdcp 

2l7 

u(r,0 )=l/7r { / /(<p) [l/2 + X ? T=i (t/ a) n cos nd cos ncp + sinnd s\nncp]dcp } 

u(r,0 )=1 In { f Q f(<p) [1/2 + En=i( r / a ) n cos n (cp - d)]dcp } ...(1) 
letc= 2^r=i( r / a ) n cos n(<^3 — 0) 

S = Y.n=iix/ay sin n(cp — 0) 

C + is = 2“ =1 {(r/a)e^- 0 )} n 

S = r/a e l ^~ e ^ / 1- r/a e l dP ~ e ); r/a < 1 and I e^ -61 ) | < 1 
S = r/a [cos {cp — 0) + i sinl'^ip — 0)] /1- r/a ( cos(<p — 0) + i sinl'jpp — 0)] 
Equating real and imaginary parts 
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C = r/a [cos (<p - 0) ± r V a 2 ] / H2 r/a) ( cos(<p - 0) + r V a 2 ] 


Using in (1) we get 


u(r,0 )=^{/; 7r /(^)[l/2 + 


27T 


-cos (cp—0)+ —7- 

1 ---^-rli 

1 “ 2 (^) cos (<p -6>) +^- 


f-27r 


(a 2 —r 2 ) 


:/(<?) 


^ ^ 271^^0 a 2 — 2 arcos((p— 9 )' 

This is the Poisson’s integral formula for a circle. 


Exterior Dirichlet Problem for a circle: 

The exterior Dirchelet’s problem is described by 

VV = 00 < 9 < 2n 

with (j){a , O) = f{0) 0 < 0 < 2n at r=a 

where /(0) is a continuous function of 9 on the surface r — a and cp must be 
bounded as r 00 . 

The solution is of the form 

Ar, 0) = z (c„ r + Z)„ r~") (A„- cos «^ + 5 ,,' sin «0) 

rc =0 

As r-> 00 ; 0(r, 9) exists finitely 

r =0 Vw 

ft 

0( r > 0)= £ r"(A, cos ” <? + sin n e ) 

n =0 


=■^ + Z r"(<3„ cos f^+b„ sin n e ) 

^ n =1 

Now by the given condition 
f(0)= Xa""(a„cosn0 + ^ n sinn0) 


...( 1 ) 


• ••(2) 
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* 2 it 


a„ = J/(^) cos (”(^)V^ 


In 


b n ~ - J/(^) Sin (”(^)V^ 

n 0 

Therefore we have from (1) 


f n\ 

a 


yr j 


— [cos n 0 cos n y/ + sin n 0 sin n ^]/(^)j y/ 


\ 2n i oo / \ w 2;r 1 

</>{ r ’Q)= — \f(v / ) d V / + -Tj T j -[cosn(^-^)]/(^)j^ 

A TC n= \ \1 r J q 


\ 2.71 1 CO / 

^• 0 ) = -f/Wi + Zr [cos(^-^)]}j^ 

o ^ n=l w 


V_ 


LetC=V - 


cos n 




And S = 


TV* - , 


sinn 




00 

T a l 


Therefore C + iS = V 

e i( ^) 

n =1 

_UJ 



$ i(y/-0) 

e 

r _ 

. ® i(^-0) 

1_ £ 
r 


Now by rationalising and comparing real parts on both sides we get 


• ••(3) 
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-cos(y/ -6)-^r 

r _r_ 

2 

1-—cos (y/-0) + ^r 
r f 


-cos(i '//-0)-Q^ 

1 ^ 71 \ V 

(/)(r,6) = -^f(xi/\- + --( 4 ) 

"° l-^cos^-^ 

r r 

i 2 r (f-a)f(v)dy/ 

'In 0 y 2 -2arcos(y/-6)) +cf 

INTERIOR NEUMANN PROBLEM FOR A CIRCLE 

The interior Neumann problem for a circle is defined as follows: 

To find the value of U at any point in the interior of the circle r=a such that 

= 0; 0< r < a; 0< 6 < 2n 

du du(r,0 ) 

And — = —--- = g(0) on r=a 

dn dr 

By the method of separation of variable the general solution of the given equation is 
given by 


= Z (C r + D n r~")(A„ cos ni9 + B n sin n0 ) 


At r=0 the solution u should be finite and therefore D n =o V n 


Therefore u(r,0 ) = ^ v\d n cosn0 + Jy sin n0) 


u(r, 0) = ^r + tl r"(a „ com6 + b„ sin n0 ) 

" n=l 

Therefore — = ^nr n ~ l {^ cosn0+ f) sinn0) 

dr n=l 
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Now ^ u ( r ’ 6) _ where 
dr 


a n 


Y 2/r 

= —— fgWcosn^ 
na n J 


2n 

b n = n -1 fgWsinn^ 

na n J 


Therefore 


u(r,0 ) = —2 + V-zr f g(^)[cosn^cos«# + sinn^sinft6?]c/0 

2 zrt nm n 


Home Assignment 

Exercise: By separating the variables, show that the equation V 2 v = 0 has a solution 
of the form A exp (+nx + iny) where A and n are constants 

Deduce that the function is of the form 


Z -mx ,rny\ 

A r e a sin(-1 x>0; 0<y<a 

r 

Where A r 's being constants are plane harmonic functions satisfying the conditions 

v(x, 0) - 0, v(x, a) — 0, v(x,y ) 0, as x -> oo 

Exercise: A thin rectangular homogeneous thermally conducting plane occupies the 
region 0 <y<b, 0<x<a. The edge y = 0 is held at temperature t(x)(x — a), 
where T is a constants and other edges are maintained at 'O'. The other faces are 
insulated and there is no heat source or sink inside the plate. Find the steady state 
temperature inside the plate. 

Sol. Pleases try yourself. 


PARABOLIC DIFFERENTIAL EQUATIONS 

The have equation of the form Rr + Ss + Tt + /(x, y, z, p, q) — 0 with S 2 — ART = 

0 is known as parabolic differential equation. The diffusion phenomenon such as 
conduction heat in solids and diffusion of viscous fluid flow as generated by a PDE of 
parabolic type. 
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The general equation for heat transfer is governed by the following equations 

dT 

ar = kVT 

Where — is the time derivative and V 2 = represents the derivative 

dt dx 2 dy 2 dz 2 r 

w.r.t., space. 


Heat Equation: The heat conduction equation = k 


May have numerous solutions unless a set of initial and boundary conditions are 
satisfied. The boundary conditions are mainly of three types and are briefly given 
below. 


Boundary condition I: The temperature is prescribed all over the boundary surface. 
This type of boundary condition depends on the problem under investigation. Some 
times the temperature on the boundary surface is a function of position only or is a 
function of time only or a constant. A special case includes T(r, t) = 0 on the surface 
of boundary, which is called a homogenous boundary condition. 

Boundary condition II: The flux of heat, i.e. the normal derivative of temperature 
is prescribed on the surface of boundary. This is called the Neumann condition. A 
special case includes =0 on the boundary.This homogenous boundary condition is 

also called insulated boundary condition which states that the heat flow across the 
surface is zero. 


Boundary condition III: A linear combination of the temperature and the heat flux is 
prescribed on the boundary 

i.e. K ^ +ht = G(x, t) 

on 

this type of boundary condition is called Robins condition. It means that the boundary 
surface dissipates heat by convection. By Newton’s law of cooling, we have 

= h(T - T a ) 

on 

T a is the temperature of surrounding 
Its special case may be taken as 

dT 

K— + hT = 0 
dn 
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Which is homogenous boundary condition. 

The other boundary conditions such as the heat transfer due to radiation obeying the 
fourth power temperature law and these associated with change of phase like melting, 
ablation etc. give rise to non linear boundary conditions. 

SEPARATION OF VARIABLE METHOD: 

We consider the one dimensional heat conduction equation 



dT , d 2 T 

dt dx 2 

...( 1 ) 

let 

T(x, t) = X{x) Y(t ) 

• ••( 2 ) 

be the solution of the differential equation ( 1 ) substituting from ( 2 ) into ( 1 ) we obtain 

X" 

x ~ 

1 Y' 

— — = A (separation parameter) then we have 

K Y 



dZz 2 " AX = 0 

ax z 

...(3) 


dY 

— - KAY = 0 

dt 

• ••(4) 


In Solving equations (3) and (4) three distinct cases arise. 


Case I: Let-A>0, say a 2 the solution will have the form 

X = Cie“ x +C 2 e~ ax , Y = C 3 e“ 2fet ...(5) 

Case II: let A = - a 2 , a is positive, then solution will have the form 
Which provide us X = Ci cos ax+ C 2 sin ax, and Y = C 3 e~ a2kt ...( 6 ) 

CASE III: let A =0 then we have 

X = C 1 x + C 2 ,Y = C 3 ...(7) 

Thus various possible solutions of the one dimensional heat conduction equation (1) 
are 

T (x, t) = (Ae ax +Be~ ax ) e fe “ 21 

T(x, t) — (A cos ax + B sin ax ) e ~“ z kt ...(8) 

T(x, t) = {Ax + B ) where A = Ci C 3 , B = -C 2 C 3 
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Example 2: Show that the solution of the equation 

dT _ d 2 T 
dt dx 2 

Satisfying the conditions 


...( 1 ) 


(1) T -» 0, as t -> oo 

(2) T = 0, for x = 0 and x = a for all t > 0 

(3) T = x, when t = 0 and 0 < x < a is 

T(x, t) = 2a /n'Zn=i(— l) n_1 /n ) sin (^- x) exp[-(n7r/a) 2 1] 
Solution: we know that the solution of (1) is 

[1] T(x,t ) = (Ae ax + Be~ ax ) exp (art) 

[2] T (x, t) = {A cos ax + B sin ax ) exp (crt) 

[3] T(x,t ) — Ax + B 


Clearly solutions represented by (1) and (2) does not satisfy the given conditions. 
Therefore the most feasible solution for the equation (1) can be treated (2) 

T(x,t ) = (A cosax + B sinax)e _a 1 using the boundary condition (2) we have 

0 = [A(l) + B(0) ] e~ a2t 

Or 0 = Ae~ a2t 


Or A = 0 

Also T(0,t) = 0 = (B sin aa)e~ a2t 
Since B =£ 0 and e _ “ 2f A 0 


=> sin aa = 0 

=>aa =n7r or a =n7r /a 

Hence the solution is of the form 

T(x, t) = B sin e~ a2t 

=Bsin (if x ) exp ( _ ^~ t 
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Since the heat conduction equation is linear therefore the most general solution is 
obtained by applying the principle of superposition 

i.e. T(x, t) = Y,n =1 Bn sin x) exp(- 

now using condition (3) we get 
T = x for t = 0 

=> x = = £“=i B n sin x 1 ( v 1 = 0) 

Which is a half range Fourier sine series therefore B n = 2/a / Q a xsin(^x) dx 

T t 7lU 

Let — x = z 

a 

— dx =d z 

a 

For x = 0, z= 0 


For x = a, z = n 7 r 

Therefore B n — - f™ — 7 zsinzdz 

a nn z 

_ 2a -l n+1 
n n 

rr^i r rns n 2a r , O0 (— l) n+1 . fnn \ ( n 2 n 2 \ 

Therefore T (x, t) — — S n =i —-—sin l — x) exp (—— t) 

EX: The ends A and B of a rod, 10 cm in length are kept at temperature 0°c and 100° 
c respectively until the steady state conditions prevails. Suddenly the temperature at 
the end A is increased to 20° c and the end B is decreased to 60° c. Find the 
temperature distribution in rod at time at t. 


Sol. The problem is described by 


dT _, d 2 T 

dt dx 2 


; 0 <x <10 


Subject to the conditions 


T (0, t) = 10 
r(10, t) = 100 

d 2 T 

For steady state —7 = 0 
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Which implies that T s = Ax + B 

Now for x = 0, T= 0 implies that B = 0, Therefore T s = Ax 
And for x = 0, T = 100° C, implies that A = 10 
Thus the initial steady temperature distribution in rod is 
T s (x) = lOx 

Similarly when the temperature at the ends A and B are changed to 20°c and 60°C, 
the final steady temperature in rod is 

T s (x) = 4x + 20 

Which will be attained after long time. At any instant of time the temperature 
T (x, t) in rod is given by 

T(x,t) — T t {x, t) + T s {x) 

Where T t (x, t) is the transient temperature distribution which tends to zero as 
t -» oo. Now T (x, t) satisfies the given partial differential equation. Hence its 
general solution is of the form 

T(x,t) — T t {x, t) + T s (x) 

T(x, t) = 4x + 20 + e~ K ^ 2t (B cos Ax +Cs in lx) 

For x = 0, T = 20°C, we obtain 

20 = 20 + B e“° 2t => B = 0, t>0 

For x = 10, T = 60° we get 

60 = 60 +e _/a2f C sin 10 A 

=> sin 10 A = 0 => % = —, ne I 

10 

The principle of superposition yields 

oo 

T(x, t) = 4x + 20 +^C„ exp 

/ 7=1 

using the initial condition T = 10 x, when t= 0, we obtain 




2 



- fcl 

[ rur 1 

t 

sin 

1 nn 
— x 


vlo"; 



U0 J 
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lOx = 4x + 20 + 7,C„ sin 

n =1 


^ YITC ^ 
-X 

V 10 y 


Where C n = — f (6x - 20)sin 
10 Jo 


r nn 3 
— x 

V 10 j 


dx 


-1 

T 


(-i) 


r 800 200 


nn nn 


Thus the required solution is 


T(x,t) = 4x + 20- [(-l) 

n= 1 L 


' 800 200 


/z/r 


tin 


exp 



f \ 

2 


/ \ 

k 

nn 

t 

sin 

nn 

— X 


Jo~J 



L io J 


Diffusion equation in cylindrical coordinates 

Consider a three dimensional diffusion equation 


8T 

8t 


KVT 


In cylindrical coordinates (r, 9, z) it become 
18T^ 8 2 T 

k 8t dr 2 r dr r 2 86 2 8z 2 


...( 1 ) 


We assume separation of variables in the form 

T(r,9,z ) = R(r ) 0(0)Z(z)0(t) 

Substituting this in (1), we get 

R" 1 R' 1 0" z" 1 $ ,2 

R r R r 2 0 z k <p 1 

Where -A" is a separation parameter. 

Then (/>'+ kA<f> 2 =0 ...(2) 


1 0 " 0 2 z" 2 

=--=-g (sa y> 


R" 1 R' 

- 1 - 1 Z-- 

R r R r 2 0 
The equation in Z, R and 0 becomes 
z" - r 2 z = 0 

R" 1 R' t. 2 

— +-+ U + u~ I r~ = - 

R r R K ’ 0 

Therefore 


0 " 

2 = - — = Q 2 
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0" + Q 2 0 = 0 

i f 2 Q 2 ^ 

And R" + -R'+ f + u-— R = 0 
r { r ) 

Equations (2) and (4) have particular solutions of the form 

<t> = e 

0 = (CcosQ# + .DsinQ#) 

A /JZ -JUZ 

Z=A e +B e 


...(4) 


Equation (5) is Bessel’s equation of order Q and its general solution is 


R (r)=C l J a [yla--jU> ) + C 2 y o [VW + ^> J 

Where J ( f)and y ( r ) are Bessel functions of order Q of first and second kind 

respectively. Equation (5) is singular for r = 0, the physically meaningful solution 
must be twice continuously differentiable in 0< r < a. 

Hence equation (5) has only one bounded solution 

i ,e ' R «= JnUK^-fJ 2 }- 

Finally the general solution of equation (1) is given as 

T(r ,6,z,t) =exp iy-K A g' + B g f ‘~ (C cosQ.0 + DsmQ.0) j - JU ) r 


Assignment 

EX: Find the solution of the diffusion equation 

— = KV 2 T 
dt 

Ex: A uniform rod of length 1 with thermally insulated surface is initially at 
temperature 9 = 9 0 At t=0, one end is suddenly cooled to 9 = 0 °C 
And subsequently maintained at this temperature, the other end remains thermally 
insulated. Find the temperature distribution 6(x, t). 

EX: Find the solution of the 1-D diffusion equation satisfying the following 
conditions 

(i) T is bounded as t oo 

(ii) f- =0, Vt 

OX Y—r 


0 , Vt 
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(iii) T(x, 0) = x(a — x); 0 < x < 0 

EX: Solve the boundary value problem 


du _ d z u 

~di =CC 


Subject to the conditions 

du (0,t) 


(i) 

(ii) 


dx 

du (1,0) 
dx 


= 0 


= 0 


dx 2 ' 


(iii) u{x, 0) = x 

EX: Solve the following equation 


0 < x < t 


du d 2 u 
dt dx 2 

Subject to the conditions 

(i) u(x, 0) = Ssinnnx 

(ii) u(0, t) — 0 — u(l, t), 0 < x < l, t > 0. 


EX: Find the solution of the equation 


du d 2 u 
dt dx 2 

Subject to the conditions 

(i) u(x, 0) = 3sinn7rx 

(ii) u(0, t) = 0 = u(l, t), 0 < x < l, t > 0. 

EX: Find the solution of the equation 

dv d 2 v 
dt ^dx 2 

Subject to the conditions 

(i) v — Vq sin nt where x = 0 Vt 

(ii) v = 0 x oo 


HYPERBOLIC DIFFERENTIAL EQUATIONS: 

One of the most important and typical homogenous hyperbolic differential equation is 
the wave equation of the form 


Where C is the wave speed. 
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The differential equation above is used in many branches of physics and 
engineering and is seen in many situations such as transverse vibrations in strings or 
membrane, longitudinal vibrations in a bar, propagation of sound waves, 
electromagnetic waves, sea waves, elastic waves in solid and surface waves in earth 
quakes. 

The solution of wave equations are called wave functions. 

Remark: The Maxwell’s equations of electromagnetic theory is given by 
V.E = AnP 
V.H = 0 


C dt 

-y Am 1 dt 

C C dt 

Where E is an electric field, p is electric charge density, H is the magnetic field, C is 
the current density and C is the velocity of light. 


Exercise : show that in the absence of a charge, the electric field and the magnetic 
field in the Maxwell’s equation satisfy the wave equation. 


Solution: we have 


Curl E = V x E 


-1 dH 
C dt 


Consider V x (V x E) - V: 


-1 dH 
C dt 


-1 d_ 
C dt 



This implies Vx(Vx£)= —\ 

C dr 

But Vx(Vx£) Can be expressed as 
V(V • E)- V 2 E = V(A;p)-V 2 E 

= -v 2 e 


-v 2 e 


-1 d 2 E 

C 2 dt 2 
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Or M = C 2 v 2 £ 

dt 2 

Which is a wave equation satisfied by E 


Similarly we can observe the magnetic field H satisfies the wave equation 

d 2 H 


dt 2 


= c 2 V 2 H . 


Solution of Wave equation: (Method of separation of variables) 

We have 


d 2 u 2 3 " u 


= c~ 


St 2 dx 2 

Let U(x, t) = X(x) T(x) be the solution of (1) 


d 2 U 


= XT" 


" dt 2 
Using in (l),we get 


and 


d 2 U 

~dl 


= TX" 


XT" = C 2 TX" ^ — = C 2 — = A (say) 
T X 


Where A is a separation parameter 

=>T"-AT = 0 
And C 2 X"-AX =0 


T = Ae^ +Be 

CASE I: If /l > 0 say /l = 

Therefore T = Ae kt + Be k ' 
Similarly 
C 2 X"-K 2 X= 0 

Ex -E x 

X(x) = D e c +E e c 


...( 1 ) 


•••( 2 ) 

•••( 3 ) 


•••( 4 ) 
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:.U(x,t) = (A e K> + B e Kt ) (D e c x + E e ~c x ) 

This is the required solution 

Case II: if 1 = 0 

Then T" = 0 and C 2 X"= 0 


=>T = Ct + D and X=Ax+B 
Therefore U{x,t) = ( Ax+B)(Ct + D ) 
Case III: If 1 < 0 say 1 = - K 2 

r-p tr /~< 2 y ff 

^> — = -X 2 and --= -X 2 

T X 


^T" + K 2 T = 0 and C 2 X" + X 2 X 
SoT = (A cos Xt + Z? sinXt) and 


= 0 
X = 


f K K ] 

l C C 


Therefore U(x, t) = (Acos Kt + B sin Kt)\ D cos — x + £ sin — x 

l C C 

REMARK: From the above solutions of the wave equation for 0 < x < 
Subject to the conditions 
1 /( 0 , t) = 0 ; t > 0 , 

t/(Z,t) = 0 

Using the conditions in case I 
0 = U (0,t) = (Ae K ‘ +Be K ')(D + E) 

=> D+E = 0 
Now U(l, t) =0 


U(l, t )=[Ae K * + Be K> ) 


f A, A, A 

£>e c + £e c 


0 


V 


7 


Z and t > 0 


...(5) 


Department of Mathematics, University of Kashmir, Srinagar-190006 


50 



Draft PDE Lecture Notes 


Khanday M.A. 


De c +Ee c 


=0 


2 —/ 

Dp c+E = 0 


=1 By comparing coefficients in (5) and (6) 


=>2 —/ = 0=>—7 = 0 
C C 

Either / = 0 or — = 0 
C 

Therefore solution in case (1) is not acceptable 
Now using in case II we get 
0=U (0, t) =(C t+ D)B 
^B = 0 

AndO = U{l,t) = (Al + B\Ct + D) 

(A/ + B\Ct + D ) =0 
Implies A =0 
Implies A=0=B 

Now using the conditions in case III 

U(x, t) = (AcosKt + B sin Kt) (^Dcos^x + E sin 

Now U (0,t) =0 

=^> (A cos Kt + B sin Kt)(D) = 0 

=* D =0 

Also 0 = U (l, t) =(AcosKt + Bsin Kt)(Esin —1^\ 


...( 6 ) 
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.Esin — l 
C 


= 0 


f-'l 

lc J 


= nn 


K nn 

c~T 


,U(x,t) = (AcosKt + BsmKt) [ Esin — x 

1 J 


nn 


Therefore by using superposition principle 


^nn ^ 
x 

1 


US X ^ = Y,E,™ ~r x A co i~r t+ B„ &in \ 

V 1 A \ 1 J 


^ cnn^ 


^cnn^ ^ 


V t j 


Ex: By the separation of variables, show that one dimensional wave equation 

d 2 Z _ 1 d 2 Z 

aA“ A aA 


Has solution of the form A exp(± inc ± inct) 

Where A and n are constants. Hence show that the function of the form 


Z(x, t) 


=£ A, 


cos 


f rnct^ 


+Z?, sin 


^ rnct^ 


a 


sin 


f rnx^ 


JJ 


\ a ) 


Where and ff are constants, satisfying the wave equation and the boundary 
conditions 

Z(0, t) = 0 = Z{a,t) ) t > 0 
Ex: Obtain the solution of the radio equation 


a 2 v _ d 2 v 

— EL/ 


dx l 


dt 1 


Appropriate to the case when the periodic e.m.f. V 0 cos (pt) is applied at the end x=0 
of the line. 

Exercise: A tightly stretched string with fixed end points x = 0 and x = l is 
initially in a position given 


y = y 0 sur 


^ nx^ 


V l J 


It is released from rest from this position. 
Find the displacement y(x,t ) 
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Sol. We have the wave equation 

d 2 y _ 1 d 2 y 

~D^ 2 ~~C 2 W 


Such that y(0,t) = 0 = y(l,t) 


And y(*,0) = y 0 sin 3 
And at t=0 


^ 7DC^ 


V I J 


dt 

Let y(x, t) = X(x) T(t) be the solution 
Then 


5 y =T ^_ = TX » 


Dx 1 


dt 2 


and 


d2y = x d -X=xr 


dx z 


dt 2 


We have 


TX" = \xT" 
C 2 


Therefore 

y(x,t) = (AcosXx + BsmXx\CcosXct + DsmXct) 

Now, 
y(0,t) = 0 

=> A(Ccos/lct + Z)sin/lcr) = 0 


A = 0 

and y(x,t) = (B sin Ax\C cos Act + D sin Act) 
dy(x, 0) _ Q 
dt 

And 
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( dy\ 

0= — = [(Acos Ax + B sin Ax\AcD cos Act - Ac C sin/tcf )], =0 

\dtJ ,=o 

=^> (A co s Ax + B s in Ax\AcD) = 0 

^D = 0 

Thus 

= EsmAxcosAct 
Where E — Be 
=> iisin/l/cos/lcf = 0 

=> sin/l/ = 0 

=> M=nn 

. nn 

=> A = — 
l 

, . ^ t-i • ( n7r ^ r ^ 

y(x,t) = y : /7sm —x cos —cf 

n ” W y W y 

By the given condition 

y(x,0) = y 0 sin 3 y* 

v * 7 


Therefore 



We know that 


. , 3sinx-sin3x 

sin x =- 

4 


Or 
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y o 


3sin 


^ 70^ 


V i J 


-sin 


^3 7VC^ 


V i J 


- £, sin 


^ 7UC^ 


+ E 2 sin 


f 27UC^ 


V i J 


+ £ 3 sin 


V i J 


^3 70^ 


V i J 


Now comparing the coefficients on both sides we get 

£,=^,£2=0, £3=^, £ 4 =£ 5 =...= 0 


y(x,t) = 


3^0 



[ 7DC ^ 



y 0 ■ 1 

2 3 7CX ^ 


(tt\ 

sin 

— 

cos 


let - -sin 


COS 



V l ) 


l/J 

1 4 1 

{ l ) 


l/J 


let 


This is the required solution. 

PERIODIC SOLUTION IN CYLINDERICAL COORDINATES: 

In cylindrical coordinates with u depending only on r. The one dimensional wave 
equation assume the form. 



r dr v 


du ' 
dr j 


1 d 2 u 

c Y ~dS 




Assume that 


U = F(r) e hvr 


Acts as a solution 


du , iwt 

~fr =F(r) e 

~ r^- = rF\r) e a 
dr 


dr 


f du ^ 
r — 

V drj 


= F'(r) e m + rF”(r) e n 


d u 2 77 / \ l v 

-gr = ~ w F(r)e 


Now substituting in (l)we get 
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\{nr) e M + 


fV) + m + 4 F(c) = 0 

r c~ 


Which is a form of Bessel’s equation and hence we have 

F = A J 0[ 

In complete form we can write this equation as 

F =a 

Therefore the complete solution for the periodic function is 


r wr 


' wr 2 

— 

+ By 


V C ) 

•J 0 

l c J 


r , 

r wr^ 


r wr 2 



r wr A 


r wr 2 

i 

o 

V C J 

+i y 0 

\ c y. 

+ C 2 

Jo 

v c J 

~ l y 0 

v c 


U = \ A Jc 


' wr' 


V C ; 


+ B 


y c 


' wr' 


V C J 


e 


Cauchy problem for inhomogeneous wave equation: 

The wave equation 

d 2 u 2 d 2 u 
—t = C —t- + j (x, t) 

St 2 dx 2 

Subject to the initial conditions 

s 2 \ du(x,0) 

u(x,0) = Tj(x), —-— = v(x) 


dt 


d U d u 
and TJ- C ^ 


Subject to the homogenous initial conditions 


u 2 (x,0) = 0 and 


du 2 (x, 0) 
dt 


= 0 


Integrating equation (2) over the region we get 


• ••(2) 
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11 


d u 

~dP 


-C 


2 


dx 2 


dxdt = || f(x,t)dxdt 


Using Greens theorem in plane we get 


-11 


dx 


c 


2 d U 2 
dx 


A d , a > 


dt 


d Ui 
V dt J 


r' 


U~ - ^ 


2 U Ul 


dx 


dt + r 

y dt dx j 


dxdt = || fix, t)dxdt 

R 

|| f(x,t)dxdt 


Where 2R denotes the boundary of the region R. The boundary 8R comprises of 
three segments PB, PA and AB 

dx 


Along PB, 


dt 


= —c 


And along PA, 


dx 

dt 


= c 


Using these, we have from equation (4) 

^ ^ du-, , du 


l c 

BP V 


du 0 , du , , 

dt-\ - -dx 


dt 


dx 


~ C | 

J PA V 


|c —^dt + —^dx = ^f(x,t)dxdt 


dx 


We know that for any function Z=Z(x, y) 

dz , dz , 

dz = —dx H-<rv 

dr dy 

=> | Cc/w 2 - | Q/w 2 = || f(x, t)dxdt 

BP PA R 

=> Cu 2 (p)- Cu 2 (b)~ Cu 2 (a)+ Cu 2 (p) -1| f(x, t)dxdt 

R 

Using conditions given in (3), we have 
u 2 (A) = u 2 (b) = 0 

Therefore we have 

2 Cu 2 ( P) = || fix, t)dxdt 


U x+c U~ ct 


1 0 1-0 

u 2 (^ > ) = 2^’l |/(*T)dxc/i 


0 x—cf +ct 
to 


(4) 
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Where P(xo to) is any arbitrary point. From the initial conditions associated with the 
homogenous system, we know that 


| -j^ x+ct 

Ml(*> r ) = - \n( x + ct )+ -ct)]+- | v{g)dg 


Hence the complete solution of the inhomogeneous wave equation in one dimensional 
system is given by 
u(x,t ) = M| +u 2 

Two Dimensional wave equation 

52 .. a 2 


d~u d~u _ 1 8 u 
~d^ 2 + ~df~C 2 ~dt 2 

i(x,y,t) = x(x ) F(y)r(r) 


...( 1 ) 


Let 


Be the solution of the of the above 2-D wave equation 


Now 

d 2 u 

dx* 


= YTX", ^Ji = XTY", ^4- 
dy dr 


= XYT" 


Using in (1) we get 

YTX" + XTY" = \ XYT" 

C 2 

Dividing throughout by XYT we get 

x" r j_7^ 

X Y ~ C 2 T 


This will be true when each member will be a constant 


Choosing the constant suitably we get 
d 2 X 


dx 2 


+k 2 x = 0 


And 
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So that 

^ + (k 2 +l 2 )c 2 t = 0 

dt 2 V ^ 

The solutions of these equations are 

X = (2 X cosAx + (2 2 sin Kx 

Y =C 3 c osL .v+C 4 sin L y 

And 

r=c,™{i/F+z;)«+ c«™(i/a 2 +l 2 }* 

Hence the solution of two dimensional wave equation is 

w(x, y,t = (Ci c os&c+ (^jSin kxf(J cosLy + Ly) 

[cs<4Jk 2 +L 2 }< + C 6 ™{i!k 2 +L 2 }') 

D'Alembert’s solution of one dimensional wave equation: 

Consider the IVP of Cauchy type described as 


8t 2 c Ox 2 

Subject to the initial conditions 

U(x,0) = r/(x), 


d~u id 2 u 

-- q — r ;-co < x < oo, t > 0 


du(x, 0) 
dt 


=v(x) 


..(1) 

• ••(2) 


Where the curves on which the initial data r/(x) and v(x) are prescribed on the 
x — axis. The functions ?/(x) and v(x) are assumed to be twice continuously 
differentiable. 

We know that the general solution of the wave equation is of the form 

u(x,t )= f(x + ct)+g(x-ct) ...(3) 

Where f and g are arbitrary functions 
Using the given conditions 


U (x, 0) =? (*) = /(*) + *(*) 
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^^-= w ( x )= C[f'(x)-g'(x)] 
ot 

i.e. C[f'(x)-g'(x)] = v(x) 

integrating (5),we get 

/(x)-g(x) = -jv(s)ds 
c o 

Now adding (4)and (6),we have 
f{x) = T ^Y~ + ^v{s)ds 
Also subtracting (6) from (4), we have 
g(x) 


Khanday M.A. 


r/ix) 

-f 

2 

2c 




...( 6 ) 


Substituting in (3) we get 

^7j{x + Ct) 1 
2 2c 


U (x, t) 


-t-cr 

J v(s)ds 


J 


7j(x ~ Ct ) 1 

2 2c 


:—ci 

Jv(5 , )di' 


This is known as D'Alembert’s solution of one dimensional wave equation. 

Note: If v=0 i.e. the string is released from rest, the solution takes the form 

= v(x + ct) + n (x-ct) 

2 

DUHAMELS PRINCIPLE FOR WAVE EQUATION 

STATEMENT : Let R 3 be the three dimensional Euclidean space and x = (xj, x 2 , x 3 ) 
be any point in R 3 . If v = v(x,/./!.) satis lies for fixed X the partial differential equation 
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S V 9 9 

-—r--C 2 V 2 V = 0 
dt 2 


with the conditions 


v(x,0,/l) = 0 


— (x,0,X)=F(x,A) 
dt 

Where F(x,A ) denotes a continuous function defined for x in R 3 

t 

ifu(x, i) = | v(x, t - A, A)dA 


Be any continuous function, then it satisfies 

xry - C 2 V 2 v = F(x,A) 
dt 

xefi 3 , t >0 


u(x, 0) = 0 = 


du(x,0) 


Proof: We are given that V satisfies the wave equation 

d 2 X 7 7 

L_4-C 2 V 2 v = 0 


With the conditions given in (1) 


Also for 


u(x, 0 = | v{x, t - /i, A)dA 


To be the solution of (2) where v(x, t-A, a) is one parameter family solution of (*) 
Also v(x,0,^)=0 for t=X 

Differentiating eq. (3) w. r. t., t under the integral sign and using Leibnitz rule we 
have 


= v(x,0, A) +J —(x, t - A, A)dA 


Differentiating (4) again w. r. t., t we have 


—y = y f (x,0,^) + J — j- (x/-A,A)dA 

kJ L q (_/ L 
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using (*) 

Finally we have from eq. (3) 

y 2 u = \y 2 v 
0 

=> 

=■ 4^-C V <i-rixM.A, 

Ot 


Clearly it(x, 0) = 0 and 
du(x,0) _ Q 
dt 


The function v(x, t, A) is called the pulse function or the force function. 

Exercise: A rectangular membrane with fastened edge makes transverse vibrations. 
Explain how a formal series solution can be obtained. 


Sol: The given equation in 2-D is given by 


d 2 u 

dt 2 


= c L 


d 2 u d 2 u\ 

■E? + asr~ x - a 


0 < y < b 


Subject to the boundary conditions 

it(0, y, t) = u(x, 0, t) — u(x, b, t) — u(a, y,t) — 0 
And initial conditions 


...( 1 ) 


Let u(x,y,t) — A(x)T(y)T(t) 


u(x,y, 0) = f{x, y) 
du(x, y, 0) 


dt 


= g(x,y ) 


be the solution of (1) 
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d 2 u " d 2 u " d 2 u n 

Then —= X YT,-^ = XY T, — XYT 

dx 2 dy 2 dz 2 


Substituting theses in (1) we get 


1 T" X" Y" . 

-t = t + t = - a ^ 


Then r"+A 2 c 2 T = 0 


=> T(t) = E cos (Act) + F sin(Act) 

... (2) 

and £+■!=-* 


X" Y" 

^ — + A 2 = -y = E 2 (say) 


X" o o Y" o 

Then — + A 2 - p 2 and — = -p 2 

X Y 


=> X" + (A 2 - p 2 )X = 0 and Y" + p 2 Y = 0 


=> A = A cos^A 2 — p 2 xj + B sin(^/A 2 — p 2 x) 

...(3) 

and Y — C cos (py) + D sin(py) 

...(4) 


put A = r, yjA 2 — p 2 — p, p — q in (2), (3) and (4) we get 

X{x) — A cos px + B sin px 
Y (y) = C cos qy + D sin qy 
T(t) — E cos rt + F sin rt 

Thus the solution is given by 

it(x, y, t) = {A cos px + B sin px) (C cos qy + D sin qy) ( E cos rt + F sin rt) 
Now using the boundary conditions it(0,y, t) = 0, we get A = 0. 

Also it(x, 0, t) = 0, => C — 0 

And u{a,y,t) = 0 => sin(pa) = 0 

mn 

=> pa = mn => p =- 

Also it(x, b,t) — 0 => sin (qb) — 0 
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=> qb = nn => q = — 
b 

Now using the principle of superposition we get, 

u{x,y,t ) = Em=iZ”=iUmn cos(rct) + S mn sin(rct)] sin^x [sin^y] 
(A) 

where r 2 — p 2 + q 2 — n 2 ^ 


The initial condition [using in (A) ] 


w(x, y, 0 ) = / (x, y) 


which implies 


/Oriy) = £m = l £n=l Smrt [ sin V*] [ Slll T y ] 

And also — = ,g(x,y) 

g(x,y) = cr£“ = 1 2 “ = 1 B mn [sin^x] [sin^y] 


... (B) 


where 


Bmn abcr 


u u 

f f /mn \ /Tin \ 

I g(x, y) sin — x J sin y— yj dxdy 


Hence (A), (B) and (C) give the required solution. 


Exercise: Solve the IVP described by —7 — c 2 —7 = e x , given that it(x, 0) = 

dt z dy L 

- du(x, 0) 9 

5 ’—= x 

Sol: Please try self. 
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